AO-AO»0  994 

UNCLASSIFIED 

I OF  I 

AO 

A060994 


TEXAS  TECH  UNIV  LUBBOCK  DEPT  OF  ELECTRICAL  ENOInEERINB  F/0  9/3 
THREE  0RAPHICAL  TESTS  FOR  THE  STABILITY  OF  MULTIDIMENSIONAL  DIO~ETC(U) 
1977  R OECARLOf  R SAEkS»  J MURRAY  AFOSR-74-2631 

AFOSR-TR-78-0232  NL 


4-78 


i 

• 

1 ’ 

■ 

I-.-  ' 

C •> 

0 

• 

ODC  FILE  COEl^  ADA050994 


-i 


1 

L 


security  classification  of  this  page  ftWiMi  ItiilH  Enlrtrin 


p^^.^DeCarlo,  R./Saeks 


ie- 


REPORT  DOCUMENTATION  PAGE 


HEAD  INSTRUCTIONS 
UEFOKE  COMI'UETINC.  FORM 


'R-7  8— 2 3 


2.  GOVT  ACCESSION 


TITI  P, 


Three  Graphical  Tests  for  the  Stability 
Multidimensional  Digital  Filters# 


Tjp^UTHOgXji  - 


EClP|CNT*S  CATALOG  NUMBER 


Type  OF  report  a period  covered 

fnterim 


6.  pEimmiiNu  otc;.  hlmct 


ilOD  COVEI 

sz, 


CONTRACT  OR  GRANT  NUM 


t/^AFOSR-74-2631  j 


9.  PERFORMING  organization  NAME  AND  ADDRESS 


tr 


Texas  Tech  University 
^-’Dept.  of  EE 

Lubbock,  Texas  79409 


/ 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  « WORK  UNIT  NUMBERS 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Office  of  Scientific  Research/NM 
Bolling  AFB,  DC  20332 


U.  MONITORING  AGEHf>rN|iAME  ft  AHnpPcc 


"'ng  Office) 


IS.  SECURITY  CLASS,  fof  (flit  rtportj 

Unclassified 


IS*.  DECLASSIFICATION/ downgrading 
SCHEDULE 


16.  distribution  statement  (o!  Mt  Report) 


Approved  for  Public  Release  - Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetrect  entered  in  Block  20,  If  different  from  Report) 


18.  supplementary  NOTES 


IEEE  Intern  Symp  Circuits  Sys  p665-669  -1977 

J 


19.  KEY  WORDS  (Continue  on  roveree  tide  if  neceeaary  and  identify  by  block  number) 


Nyquist  Criterion,  Multidimensional  Digital  Filters,  Stability, 
Homotopy  Theory 


BSTRACT  (Continue  on  reveree  aide  if  neceaaary  and  Identify  by  block  number) 


This  paper  discusses  three  graphical  tests  for  determining  the  stability  of 
multidimensional  digital  filters  characterized  by  an  appropriate  transfer 
function  iryseveral  complex  variables.  Each  test  is  carried  out  as  a finite 
number  of  "^yguisf^plots  in  the  complex  plane. 


\ \ ■' 
W; 


DD  1 JAN  73  1473  edition  of  ) NOV  65  IS  OBSOLETE 


iiNiri.A.qsiFiED 


A5 


SECURITY  CLASSIFICATION  OF  THIS  PACE  fWh«n  Dote  Entered) 


-1 

J 


Al-'usin®,  7 8-  0232 


Part  I — All  Papers  Except  Sessions  SI  & S3 


1977  IEEE 
INTERNATIONAL 

Symposium 

ON 

CIRCUITS  AND 
SYSTEMS 


•:  ^ 'fp.  , ; / o 

Ij  j ■ .)  S 


1 


THREE  GRAPHICAL  TESTS  FOR  THE  STABILITY  OF 
HULTIOIHENSIONAL  DIGITAL  FILTERS^ 

R.  DeCarlo,  R.  Saeks,  and  J.  Murray 


Texas  Tech  University  at  Lubbock,  Texas 


, Wilts  Section  if 

D-II  Section  □ 
bNANNCOUCi-Q  □ 

JUSTIftC^TiaN 


BY 


®s_] 

' sPtClAL' 


ftM 


ii"  - 

i .! 


Abstract 

This  paper  discusses  three  graphical  tests 
for  determining  the  stability  of  multidi- 
mensional digital  filters  characterized  by 
an  appropriate  transfer  function  in  several 
complex  variables.  Each  test  is  carried 
out  as  a finite  number  of  "Nyguist”  plots 
in  the  complex  plane. 

Introduction 

Recently  two  of  the  authors  constructed  an 
algebraic  topological  proof  of  the  Nyguist 
Criterion  (2)  (3) . The  value  of  this 
rather  sophisticated  approach  has  been 
harvested  in  generalizations  to  systems 
characterised  by  transfer  functions  in 
several  complex  variables  (1)  U)  (3),  in 
particular  multidimensional  digital  filters. 
Specifically  the  paper  Illustrates  three 
graphical  tests,  similar  to  the  classical 
Nyguist  test,  carried  out  in  the  complex 
plane,  which  determine  the  stability  of  a 
multidimensional  digital  filter  with  a 

transfer  function  H(z,,...,z_)  > B(z, ,..., 

in  X 

z^)/A(z^, . . . ,z^)  where  z^  are  complex  vari- 
ables and  A and  B are  relatively  prime  poly- 
nomials. The  purpose  of  the  paper  is  to 
consider  these  three  tests  as  applied  to 
two  different  examples. 

Background  and  Main  Theorems 

Basic  to  the  theory  is  the  2n  (r«al)  dimen- 
sional polydisc  (8)  which  is  the  0“  analog 
of  the  unit  disc  of  0.  Mathematically  the 

polydisc  p"  is 

P"  - {(Zj,...,z^)in  0"|  |zj^|<  1,  i-l,...,n} 

There  are  four  separate  notions  of  boundary 
of  the  polydisc  (1)(8).  First  is  the  usual 
topological  boundary 

b"  - 1,  i-l,...,n 

and  |Z|^|  > 1 for  some  k) 
Second  is  the  distinguished  boundary 
t"  - {(zj^,...,z^)in  - 1,  i-l,...n} 

T**  serves  as  the  multidimensional  analog  of 
the  unit  circle.  In  particular  the  freguency 


response  (7)  of  a digital  filter  is  the 
evaluation  of  its  transfer  function  over 

t". 

Thirdly  we  have 

m"  - {(zi,...,z„)in  0'*||zj|-l k-l,k+l, 

. . . ,n/  ^ ^ 

where  k ranges  from  1 through  n.  This  is 
a boundary  notion  in  the  sense  that  n-1 
Coordinates  take  on  extremal  values.  Fin- 
ally, the  last  notion  of  "boundary”  is 

H * { (Zj^, . . . ,S|^) is  0 ||z^|"l,  l>l,...,k-l; 
l*kll  *i  “ i”k+lf -...n} 

where  again  k varies  from  1 to  n.  The  im- 
portance of  this  concept  was  first  noted 
by  Huang  (5) . Later  it  was  generalized  in 
(7). 

With  these  notions  of  boundary  one  may 
prove  the  following  Theorem.  The  proof, 
however,  of  the  following  eguivalence  is 
found  in  the  references  (1) (2) (5) (7) (9) . 
Theorem  1:  Let  a causal  multidemsional 
digital  filter  be  characterized  by  a 
rational  transfer  function  in  several  com- 
plex variables.  Assume  the  numerator  and 
denominator  polynomials  are  relatively 
prime.  Then  the  following  are  eguivalent 
stability  conditions: 

(i)  the  pole  set  of  the  transfer 
function  has  a null  Intersection  with  p" 

(li)  the  pole  set  of  the  transfer  * 
function  has  a null  intersection  with  b" 
(iii)  the  pole  set  of  the  transfer 
function  has  a null  intersection  with  H** 

(iv)  the  pole  set  of  the  transfer 
function  has  a null  intersection  with  M^. 

The  trouble  with  these  conditions  is  that 
the  actual  test  is  carried  out  in  a higher 
dimensional  space.  For  example  P**  is  2n 
dimensional,  B is  (2n-l)  dimensional, 

while  h"  and  M**  are  both  (n-fl)  dimensional. 
Intuitively,  the  egul valences  of  this  theo- 
rem follow  because  the  pole  set  of  a ration- 
al function  in  several  complex  variables 
is  an  infinite  continuum  which  must  inter- 
sect the  different  boundaries  of  the  poly- 
disc if  it  intersects  the  polydisc  at  all. 
Nith  the  Intuition  gained  in  (2) (3)  (7)  the 
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authors  were  able  to  simplify  these  results 
to  graphical  tests  in  the  complex  plane. 

Vhe  following  three  theorems  are  the  fruit 
of  this  endeavor.  Before  stating  these 
theorems,  one  final  definition  is  in  order. 
The  "Nygulst  plot”  of  a polynomial,  f('), 
in  one  complex. variable  is  defined  to  be 
the  image  of  T'*’  under  the  map  f(')  where 

is  the  unit  circle  of  the  complex  plane. 
With  this  in  hand,  we  have  the  following 
three  tests.  Again  the  proofs  can  be  found 
in  the  references  (1) (2) (3) (9) . 

Theorwn  2;  Let  A be  the  denominator  poly- 
nomial  of  a multidimensional  digital  filter 
as  characterized  in  Theorem  1.  A is  a 
polynomial  mapping  0 to  (f.  Then  A has  no 

zeros  in  (i.e.  the  filter  is  stable)*  if 
and  only  if 

(i)  A has  no  zeros  on  t”,  and 
(ii)  The  Nyquist  plots  for  the  single 
variable  functions  A(l, . . . ,l,Zj^,l, . . . ,1) 
k « l,...,n  do  not  encircle  zero. 

Theorem  3;  Let  A be  as  above.  Then  A has 

no  zeros  in  P**  (i.e.  the  filter  is  stable)* 
if  and  only  if 

(1)  A has  no  zeros  on  t",  and 
(ii)  The  Nyquist  plots  for  the  single 
variable  functions  A(l, . . . ,l,Zj^,0, . . . ,0) 
k » l,...,n,  do  not  encircle  zero. 

Theorem  48  Let  A be  as  above.  Then  A has 

no  zeros  in  P**  (i.e.  the  filter  is  stable)* 
if  and  only  if 

(1)  A has  no  zeros  on  t",  and 
(ii)  The  Nyquist  plot  for  the  single 
variable  function  A(z,z,...,z)  does  not  en- 
circle zeroL 

Each  of  these  tests  has  essentially  the 
same  two  parts.  First  one  performs  the 
appropriate  encirclement  test(s);  if  zero 
is  not  encircled,  one  then  proceeds  to  check 
the  image  of  the  distinguished  boundary.' 

This  order  of  testing  (encirclement  first, 
then  frequency  response)  seems  in  most  cases 
to  be  preferable  to  the  reverse  order,  since 
much  less  computation  is  involved  in  the 
encirclement  tests;  however,  in  cases  where 
the  frequency  response  is  Icnown  a priori, 
or  must  be  plotted  in  any  case,  the  order 
is  immaterial. 

It  might  appear  that  the  third  test  (Theorem 
4)  is  the  best,  since  it  Involves  only  one 
encirclement  test;  however,  in  many  cases 
the  relative  complexity  of  the  polynomial 
A(z,z,...,z)  will  more  than  offset  this 
advantage.  Similarly,  in  many  cases.  Theo- 
rem 3 may  be  much  easier  to  apply  than  Theo- 
rem 2.  (This  is  Illustrated  in  the  first 
example) . Theorems  2 and  4,  however,  do  have 
two  advantages.  The  first  is  mainly  philo- 
sophical; these  Theorems  give  a test  for 
stability  purely  in  terms  of  the  frequency 
response  of  the  function  A,  which  corres- 
ponds closely  with  the  idea  of  the  Nyquist 
criterion  in  one  variable.  The  second  ad- 
vantage is  that  by  filling  in  the  interior 
of  the  encirclement  plot(s)  and  taking  this 
*See  note  2. 


region  together  with  the  image  of  the  dis- 
tinguished boundary,  one  obtains  the  image 
of  the  entire  polydisc,  from  which  one  can 
get  an  accurate  idea  of  stability  margins. 
(The  point  here  is  that  we  have  found  the 
image  of  a 2n-dlmensional  set — the  poly- 
disc— by  plotting  an  n-dimensional  set  and 
a 1-dimensional  set) . 

EXAMPLES 

In  this  section  we  apply  each  of  the  above 
tests  to  two  examples. 

Example  1: 


A(z^,Z2)  - 5/4  zl  zl 
+ 3 z,  z?  + 3 z?  J 


+ •,  ZjZ2 


*2  * *1  ■ *2 


+ 3 


In  this  case,  we  have 

A(z,l)  - 1/2  z^  + 13/4  z^  ♦ 9/2  z-2 

A(l,z)  - 1/2  z^  + 13/4  z^  ♦ 9/2  z - 2. 

These  polynomials  are  identical;  the  image 
of  the  unit  circle  Iseing  given  in  Fig.  1(a).* 
Since  this  curve  encircles  the  origin,  we 
deduce  Immediately  that  the  filter  is  un- 
stable; for  purposes  of  Illustration,  we 
will  carry  out  the  other  tests. 


A(z,0) 


-z*  - 2z  + 1 


A(l,z)  is  as  before  (Fig.  1(a));  A(z,0)  is 

plotted  (for  z •>  e^^)  in  Fig.  1(b).  Again 
either  plot  suffices  to  verify  instability, 
and  clearly  A(z,0)  gives  the  simpler  test. 
To  apply  the  third  test,  we  calculate 


A(z,z)  - 9/4  z’  + 6 Z-*  + z‘  - 4z  + 1 

and  the  image  of  the  unit  circle  under  this 
mapping  is  plotted  in  Fig.  1(c).  The 
relative  complexity  is  apparent;  however, 
it  again  verifies  instability.  Finally, 
we  plot  the  image  of  the  distinguished 
boundary  in  Fig.  1(d);  it  can  be  seen  that 
it  does  not  include  the  origin,  although 
it  dees  in  some  sense  "encircle”  it. 

The  second  example  shows  that  this  last 
kind  of  "encirclement"  is  irrelevant; 
nothing  can  be  deduced  from  it. 

Example  2; 

A(Zj,Z2)  - (Zj  + 2)^  (Z2  + 2)^ 

As  before,  the  plots  for  A(z,l)  and  A(l,z) 
are  identical 

A(z,l)  » 27(z  + 2)? 

A(l,z)  - 27(z  + 2)^; 

This  plot  is  given  in  Fig.  2(a);  it  does 
not  encircle  0. 

*(This  Illustrates  the  obvious  fact  that 
if  the  polynomial  is  symmetric  in  s.,..., 
Zj^,  then  the  n plots  in  Theorem  2 in  fact 

reduce  to  1 plot — usually  simpler  than  the 
plot  in  Theorem  4.  Such  symmetry  is  quite 
common) . 


- 4z  + 1 


A(Zj^,Z2) 


In  this  case,  the  plot  tor  A(z,0)>>8 (z-i-2) ^ 
differs  from  the  previous  plot  only  by  a 
scale  factor;  do  not  draw  it  separately. 
Finally,  the  plot  for  A(z,s)-(z  -f  2)^  is 
given  in  Fig.  2(b);  again  it  does  not  en- 
circle the  origin. 

Thus,  in  order  to  determine  stability  in 
this  case,  it  is  necessary  to  plot  the 
image  of  the  distinguished  boundary;  this 
is  done  in  Fig.  2(c).  Since  this  image 
does  not  contain  the  origin  (although  it 
does  surround  it) , we  conclude  that  the 
filter  is  stable.  This  of  course  is  ob- 
vious analytically;  the  present  example 
is  merely  to  Illustrate  the  tests. 

Note  It  Because  of  the  magnitudes  of  the 
numbers  involved.  Fig.  2(a)  - 2(c)  are  not 
drawn  to  scale. 

Note  2;  Nonvanishing  of  the  denominator 
polynomial  implies  stability;  the  converse, 
however,  is  not  quite  true  unless  we  assume 
that  the  transfer  function  has  no  Indeter- 
minacies  on  the  distinguished  boundary 
(see(lO)]. 
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